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ABSTRACT 

We measure proper motions with the Hubble Space Telescope for 16 extreme radial velocity stars, 
mostly unbound B stars in the Milky Way halo. Twelve of these stars have proper motions statistically 
consistent with zero, and thus have radial trajectories statistically consistent with a Galactic center 
“hypervelocity star” origin. The trajectory of HE 0437-5439 is consistent with both Milky Way and 
Large Magellanic Cloud origins. A Galactic center origin is excluded at 3cr confidence for two of 
the lowest radial velocity stars in our sample, however. These two stars are probable disk runaways 
and provide evidence for ~500 km s“® ejections from the disk. We also measure a significant proper 
motion for the unbound sdO star US 708. Its 1,000 km s“® motion is in some tension with proposed 
supernova ejection models, but can be explained if US 708 was ejected from the stellar halo. In the 
future, we expect Gaia will better constrain the origin of these remarkable unbound stars. 

Subject headings: stars: early-type — kinematics and dynamics — Galaxy: halo — Galaxy: kine¬ 
matics and dynamics — Magellanic Clouds 


1. INTRODUCTION 

The Milky Way is a gravitationally bound system 
of a few 10^® stars, but it also hosts some unbound 
stars. Where these unbound stars come from is an open 
question. It is very difficult to explain unbound main- 
sequence staxs_by_simernova explosions in close binary 
systems (lBlaaiiwlll961h or by dyn amical encounters be¬ 
tween binaries (iPoveda et al.lll9CT . the mechanisms that 
explain disk “runaway” stars. Simulations show that 
stars cannot be launched at speeds greater than their bi¬ 
nary orbital velocity, and this velocity cannot exceed the 
escape velocity from a star’s surface, since that would re¬ 
quire the stars to be orbiting inside one another (|Leonardl 
II99II1 . The escape velocity from the surface of main- 
sequence stars is comparable t o the 500 - 600 km s“® 
local Galactic escape velocity (|Piffi et al.l 120141 1. Only 
dynamical interactions with an object more massive and 
compact than a star c an ea sily explain unbound main- 
sequence stars. iHillsI (|1988ll predicted that three-body 
exchange interactions between stars and the central mas¬ 
sive black hole (MBH) will inevitably unbind stars from 
the Galaxy, and he called such objects “hypervelocity 
stars” (HVSs). We therefore refer to unbound main- 

sequence stars in this paper as HVSs^ _ 

The first HVS was discovered bv iBrown et al.l (|2005ll . 
and a couple dozen B-type HVSs are now known in 
the M ilk y Way ha M (lEdeimann et al.l l2005l: IBrown et al.l 
2006a [0.l2007allbl: iHeber et al. I2008I [Brown et al.l 120091 
20121 ' l2014HZheng et al.ll2014ll . These stars have radial 
velocities up to -1-700 km s“®; no star moving towards us 

wbrown@cfa. harvard.edu;jayander®stsci.edu;ognedin®umich.edu 
^ Based on observations with the NASA/ESA Hubble Space 
Telescope obtained at the Space Telescope Science Institute, 
and from the data archive at STScI, which is operated by the 
Association of Universities for Research in Astronomy, Inc., 
under NASA contract NAS5-26555. 


with a comparable radial velocity has ever been observed. 
To date, detailed spectroscopic analyses of the B-type 
HVSs find that th ey are main-sequence B stars at 10 - 


100 ki 

DC distances (lEdeimann et al.ll2005l:IPrzvbilla et al. 

2008a 

Ibl: iLoDez-Morales & BonanosI 120081 IBrown et al. 

2012a 

.1201311. The short lifetimes of B stars require that 

these unbound stars were ejected from a region with re- 


cent star formation, such as the Galactic disk or the 
Galactic center. 

Full three-dimensional trajectories for the HVSs are a 
crucial test of their origin. Measuring the trajectories for 
the HVSs, however, requires accurate distances and ab¬ 
solute proper motions. Distance estimates to individual 
stars are determined by comparing either spectroscopic 
Teff and log g or broadband colors to stellar evolutio n 
tracks. Typical precision is 15% (jBrown et al.l 1201^ : 
the accuracy depends on the choice of stellar evolution 
tracks (i.e. assumptions about metallicity). Proper mo¬ 
tions are equally difficult to measure, because HVSs are 
distant and they should be on radial trajectories. Ex¬ 
pected HVS proper motions are typically < 1 mas yr“® 
and thus cannot be measured with ground based sur¬ 
veys. The only unbound star with a statistically signifi¬ 
cant proper motion measurement to date is the runaway 
B star HD 271791, a star that was ejected in the direc¬ 
tion o f Galactic rotation from the outer disk (jHeber et al.l 
l2008fl . HIP 60350 is another example of a high velocity B 
star with a significant proper motion mea surement, but 
it is u nbound at less than Icr significance (|Irrgang et al.l 
|20I3). Here, we present Hubble Space Telescope (HST) 
proper motion measurements for a sample of 16 stars 
with extreme radial velocities, 12 of which are unbound 
to the Milky Way in radial velocity alone. 

Proper motions are important because they can dis¬ 
criminate between Galactic center and Galactic disk ori¬ 
gins. The star US 708 (hereafter HVS 2), for example, is 



















































2 


Brown, Anderson, Gnedin, et al. 


a helium-rich sd O star with a -1-70 8 km s“^ heliocentric 
radial velocity (|Hirsch et alJ 120051) . sdO stars are the 
remnants of low mass stars; helium rich sdO stars are 
possibly the merger product of two helium white dwarfs 
l|Hebeil2009li . IfHVS 2 were ejected by the central MBH, 
its velocity vector will point away from the Galactic cen¬ 
ter. If HVS 2 were ejected by a Type la supern ova ex¬ 
plosion (lJustham et al.l[200^ (Wang fc Hanll200^ , its ve¬ 
locity vector can point from anywhere in the Milky Way. 

Generally speaking, lower velocity stars are more likely 
to be disk runaways . The maximum ejection velocity 
in the iBlaanwl (jlQBlh supernova mechanism is the sum 
of the supernova kick velocity and the orbital veloc¬ 
ity of the progenitor binary, or about 250 km s~^ fo r 
an ejected main-sequence B star (jPortegies Zwartir200ClD . 
Velocities up to 500 km s“^ may be possible for B stars 
ejected by extremely asy mmetric cor e-collapse super¬ 
novae in contact binaries (lTaurisl[2015D . The maximum 
ejection velocity in the dynamical mechanism depends 
on the most massive star in the encounte r, and can also 
reach 500 km s~^ in ex treme scenarios (jLeonar I IWH: 
iGvaramadze et al'1l2009li . Simulations suggest that 99% 
of dynamical ejection s occur at velocities <200 km s“^ 
(jPerets fc Subd[20T^ . We test for the existence of ^500 
km s“^ runaway ejections by including four bound veloc¬ 
ity outliers in our sample. 

The star HE 0437-5439 (hereafter HVS 3), an unbound 
9 Mq main-sequence B star located 16 deg from the 
Lar ge Magellanic Gloud ( LMG), presents another puz¬ 
zle (jEdelmann et al.ll2005li . The lifetime of this star is 
many times shorter than its flight time from the Milky 
Way, suggesting an LMG origin. An LMG origin re¬ 
quires a 1,000 km s“^ eje ction, however, and thus an un¬ 
seen MBH in the LMG (|Gualandris fc Portegies ZwartI 
l2007f) . A Milky Way origin, on the other hand, requires 
that HVS 3 be a blue straggler to account for its life¬ 
time. In other words, the progenitor must have been a 
binary system ejected at > 800 km s“^ that subsequently 
merged as it traveled away from the Milky W ay. A bi¬ 
nary MBH could eject stellar binaries as HVSs (iLu et al.l 
UMl), or else a sin gle M BH could eject binaries by triple 
disruption (lPeretsll2009li . All of these models have very 
low ejection rates. The LMG and Galactic center ori¬ 
gins for HVS 3 differ by 1.5 mas yr ~^ in proper motion . 
Using two epochs of HST imaging, iBrown et al.l (j201flf j 
found that HVS 3 is moving away from the Milky Way. 
lirrgang et al.l (|2013D argue that both Milky Way and 
LMG origins, given the systematic uncertainties, are con¬ 
sistent with the measurements. A third epoch of imaging 
is needed. 

In principle, HST proper-motion measurements have 
sufficient precision to determine the origin of our un¬ 
bound stars. Using the best data reduction and measure¬ 
ment techniques, it is possible to achieve 0.0 1 pixel astro- 
metr ic precision on well-exposed stars (e.g. iBellini et al.l 
1201411 . The Advanced Camera for Surveys (ACS) Wide- 
Field Channel has 50 mas pixels, and the Wide Field 
Camera 3 (WFC3) Ultraviolet-Visible (UVIS) Channel 
has 40 mas pixels. For a pair of measurements separated 
by a 3 year time baseline, we can thus expect 0.2 mas 
yr“^ precision. We must rely on background galaxies to 
define our absolute reference frame, however, and both 
the number and spatial distribution of useful background 


galaxies introduce a systematic uncertainty to the mea¬ 
surements. In some cases, a sequence of short and long 
exposures is used to link faint background galaxies to 
bright HVSs, adding additional uncertainty. 

For the sake of clarity, we quote proper motion uncer¬ 
tainties that are the sum in quadrature of the statisti¬ 
cal (stellar) and systematic (background galaxy) proper 
motion uncertainties throughout this paper. Our average 
total uncertainty is ±0.8 mas yr“^, a 6-fold improvement 
over existing proper motion measurements. 

In Section 2 we define the sample of 12 unbound and 
four bound stars, and present new spectroscopy and stel¬ 
lar atmosphere fits for the four bound stars. In Sec¬ 
tion 3 we describe the HST imaging, image reduction, 
and proper motion measurements. Twelve stars have 
proper motions statistically consistent with zero, and 
thus largely radial trajectories. In Section 4 we evaluate 
the constraints on a Galactic center origin, accounting 
simultaneously for radial velocity, distance, and proper 
motion errors. Thirteen stars have trajectories consistent 
with a Galactic center origin given the measurement er¬ 
rors, while three stars are inconsistent with a Galactic 
center origin at >3cr confidence. In Section 5 we discuss 
the objects with significant proper motions, the runaway 
stars HVS 2, B 711, and B 733. We also update our 
constraints on HVS 3, consistent with either a LMG or 
Milky Way origin. We conclude in Section 6. 


2. THE SAMPLE 

We select our sample of 16 stars for their extreme ra¬ 
dial velocities. We targeted all 12 unbound HVSs known 
prior to 2008 March, plu s four bound velocity outliers 
from the HVS Survey of IBrown et al.l (l2007blj . Table 
[T] lists the coordinates and other observed properties of 
the stars. We refer to the 12 unbound st ars as HVS 1 
(IBrown et al.l 120051) . US 708 = HVS 2 (iHirsch et al.l 
I2005D . HE 0437 - 5439 = HVS 3 (lEdelmann et al.ll20n5D . 
HVS 4 - HVS 5 (IBrown et al.ll2006aD , HV S 6 - HVS 7 
(iBrow-n et al.l l2006b[ ). HVS 8 - HVS 10 (IBrown et al.l 
l2007bD . and HVS 12 - HVS 13 (IBrown et 311120091) . We 
refer to the four bound velocity outliers as B 434, B 485, 
B 711, and B 733, c orresponding to their target numbers 
in the HVS Survey (|Brown et al.ll2007E) . 

Detailed analyses of the unbound stars indicate that, 
except for HVS 2, they are main-sequence B stars at 
50 - 100 kpc distances. This conclusion is based on 
stellar atm osphere fits to high resolution echelle spectr a 
of HVS 3 (iPrzvbilla et al.l l2008at iBonanos et al.ll2008D. 
HVS 5 (|Brown et al.l 2M2aD. HVS 7 (IPrzvbi l la et alJ 
l2008bD . and HVS 8 (iLopez-Morales fc BonanosI I2008D . 
and fits to moderate resolution, high signal-to-nois e 
(S/N) spectra of the other HVSs (jBrown et al.ll20l3) . 
Here we describe spectroscopy of the four bound velocity 
outliers. 


2.1. New Spectroscopy 

We obtained 6.5m MMT Blue Channel spectroscopy of 
the four bound velocity outliers in four different observ¬ 
ing runs in 2006 February, 2008 February, 2010 March, 
and 2014 March. We used the 832 1 mm“^ grating in 
2nd order with a 1 arcsec slit, providing 1.0 A spectral 
resolution over the range 3600 - 4500 A. In addition, we 
observed B 485 a second time with a 0.75 arcsec slit to 
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Fig. 1. — Spectra of the bound velocity outliers, continuum- 
normalized and shifted to rest-frame (in black), compared to their 
best-fitting stellar atmosphere model (in red). The wavelength sep¬ 
aration between the pair of Ca ii A3933 lines (marked with arrows), 
the lefthand one due to local interstellar medium absorption, vis¬ 
ibly shows each star’s large radial velocity. We use the hydrogen 
Balmer lines to measure Teff, log and usini. 

test higher 0.75 A spectral resolution. We paired every 
observation with a He-Ne-Ar comparison lamp exposure 
for wavelength calibration. We chose exposure times to 
achieve S/N=50 - 100 per resolution element in the con¬ 
tinuum. 

We process and extract the ID spectra using IRAFQ 
We measure radi al velocities using the cross-correlation 
package RVSAO (|Kurtz fc MinkI flQQSll . We then mea¬ 
sure effective temperature Tea, surface gravity log 5 , and 
projected rotati on v sin i usiiig stella r atmosphere models 
as described in iBrown et al.l 1201411 . Figure [T] compares 
the best-fit stellar atmosphere models to the observed 
spectra. All four bound velocity outliers are late B-type 
stars, and B 733 has significant projected rotational ve¬ 
locity usinf=240±30 km s“^. 

2 . 2 . Stellar Parameters 

Our sample contains 15 B-type stars plus the sdO star 
HVS 2. Figure [5] displays the adopted effective temper¬ 
ature Teff and surface gravity logg for the 15 B-type 
stars in our sample, plotted in comp arison to Padova so¬ 
lar metallicity main-sequence tracks (lOirardi et al1l2004 
iMarigo et al.ll200§ : iBressan et al.ll2012l l. HVS 2 does not 

appear because it is a reff=45,560 K sdO star. The 
clumping of stars in Figure [5] reflects the HVS Survey 
target selection of stars with the colors of 3 Mq main 
sequence stars. 

Notably, all four stars with high resolution echelle spec¬ 
troscopy are fast rotators. Another five B-type stars have 
statistically significant usinf > 150 km s“^ on the basis 
of our moderate resolution MMT spectroscopy. Fast ro¬ 
tation is the unambiguous signature of a main-sequence 
B star. For reference, the mean v sin i of late B-type stars 

^ IRAF is distributed by the National Optical Astronomy Ob¬ 
servatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 


Fig. 2. — Effective temperature Teff and surface gravity logg 
compared with Padova solar metallicity main-sequence tracks (red 
lines) for our sample of stars. Stars are labeled by their HVS/B 
identifier number. 

is 150 km s“^ (|Abt et al.l l 2002 t iHuang fc GiesI I2006D . 
Evolved horizontal branch stars of the same temperature 
and surface gravity, on t he other ha nd, are slow rotators 
with usinz < 7 km s“^ (|Behiil200^ . Thus these B-type 
stars are main-sequence stars. 

We compare the measured stellar atmosphere parame¬ 
ters to Padova main-s equence stellar evolution tracks to 
estimate luminosities (iGirardi et akl 12004 IMarigo et al.l 
I2OO8I : IBressan et al.l I 2012 D . Propagating the measure¬ 
ment uncertainties through the stellar evolution tracks 
implies that our luminosity estimates are precise to 30% 
in luminosity, or 15% in distance. The precision is rel¬ 
atively poor because surface gravity is our primary con¬ 
straint on evolutionary status, and the luminosities of B 
stars increase with age. 

The choice of stellar evolution tracks is a source of sys¬ 
tematic uncertainty. With the exceptio n of HVS 3 and 
HVS 7 , which have solar iron abundance (jPrzvbilla et al.l 
l2008airbl l. the metallicity of the HVSs is unconstrained. 
If the HVSs are systematically metal-rich or metal-poor, 
our distanc e estimates could be systematically wrong by 
^^25% (e.g. IBressan et ai]l2012t) . Given that the HVSs 
are relatively short-lived main-sequence B stars, however, 
we consider solar metallicity a reasonable assumption. 

We calculate heliocentric distance using the de- 
reddened 5 -band point spread function apparent mag¬ 
nitudes from Sloan Digital Sky Surv ey (SDSS) Data Re¬ 
lease 10 (DRIO) (lAhn et al.l 12014. except for HVS 3 
which has a V-band measurement (iBonanos et al.ll200^ . 
The SDSS magnitudes have a typical precision of 2%; 
however in some cases the DRIO 5 -band magnitude dif¬ 
fers from previous SDSS data release values by as much 
as 10%. The SDSS photometry is thus an additional 
source of uncertainty, but a small fraction of the total 
error budget. 

Table □ summarizes the measured and derived prop¬ 
erties for the 16 stars. We adopt parameters from pub¬ 
lished echelle spectra when available; otherwise we adopt 
parameters derived from our MMT spectra. The stellar 
atmosphere parameters for the four bound velocity out- 
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liers are new. 

3. PROPER MOTION MEASUREMENTS 
3.1. HST Observations 

We imaged the 16 stars using the HST ACS and 
WFC3/UVIS instruments starting in 2006 September 
and ending in 2013 March. We obtained images in a few 
different ways; Table [5] lists the full set of observations. 

We first observed HVS 1 - HVS 5 using ACS in 2006 
(proposal ID 10824). Each star was allocated one or¬ 
bit, and was observed with a set of 4 - 7 dithered expo¬ 
sures. Exposure times were chosen to maximize counts 
on background galaxies while not saturating the HVSs. 
We obtained a second epoch of observations with ACS 
in 2009 (proposal ID 11782). Eor HVS 1 - HVS 4, we 
then obtained a third epoch of observations using WFC3 
(proposal ID 12503). We allocated two orbits of time per 
star in the third epoch, but dropped HVS 5 because its 
field contained few useful background galaxies. The time 
baseline of observations for HVS 1 - HVS 4 is thus about 
6.2 yrs, and for HVS 5 it is 3.0 yrs. 

We observed the other 12 stars exclusively with WFC3 
beginning in 2009-2010 (proposal ID 11589). Each star 
was allocated one orbit of time and observed with a set 
of dithered exposures. Eor the brightest 6 stars, how¬ 
ever, we divided the exposures into a set of short ~1 min 
and deep ~10 min exposures. Our goal was to optimally 
expose the bright HVSs and faint background galaxies 
in the short and deep exposures, respectively, and to tie 
their astrometric frames together using the intermedi¬ 
ate brightness stars available in both sets of exposures. 
In practice, the finite number of intermediate brightness 
stars limits the accuracy of this approach. We obtained 
a second epoch of observations with WFC3 in 2012-2013 
(proposal ID 12662) in the same way. The time baseline 
of the observations is 3.0 yrs for these 12 stars. 

We made an effort to use the same telescope orien¬ 
tation across all epochs of imaging in order to minimize 
the impact of changes in charge-transfer efficiency (CTE) 
and errors in the distortion solution. The effects of CTE 
increase with time as on-orbit radiation damage creates 
more and more charge traps in the CCD silicon lattice. 
The result is that star and galaxy images are trailed along 
the CCD read-out direction. If we observe in the same 
orientation in all epochs, the CTE systematic is in the 
same direction and we minimize its impact on our dif¬ 
ferential position measurements. Unfortunately, we were 
unable to use the same guide stars for HVS 2, HVS 3, 
and HVS 4 in all epochs due to changes in the HST guide 
star catalog. This issue forced a change in orientation, 
which exacerbates the astrometric effects of CTE. Table 
[2] lists the position angle (P.A.), defined as increasing 
East of due North, for each observation. 

3.2. Image Reduetion and Analysis 

We begin our data reduction procedure by down¬ 
loading flat-fielded images from the Mikulski Archive 
for Space Telescopes. Because the CTE correc¬ 
tion for WFC3/UVIS images and ACS subarrays is 
not integrated into the st andard pipeline, we run 
the lAnderson fc Bedinl (|2010ll pixel-based correction by 
hand. The CTE correction is calibrated on hot pixels and 
their charge trails, and is successful at removing the trails 


behind stars, cosmic rays, and hot pixels, and restoring 
the flux to the stellar images. 

We start our analysis with the first-epoch images. We 
use empirical models of the ACS and WFC3/UVIS point 
spread functions (PSEs) to measure positions for stars 
in the CTE-corrected frames of each first epoch expo¬ 
sure. We then correct the positions for geometr ic distor¬ 
tion using the mo del inlAnderson fc Kind (|2006ll for ACS 
and the model in I Bellini et al.l ( 201111 for WFC3/UVIS. 
Finally, we cross identify the stars and define a linear 
transformation from the distortion-corrected positions of 
each star in each exposure into the distortion-corrected 
frame of the first exposure. 

Since all of the first-epoch images were taken in a single 
orbit, we can safely use the stars to define the transfor¬ 
mations. This cannot be done for the later epochs, since 
the stars are all moving, and we need to know the motion 
of each target star in an inertial frame. Therefore, we 
use the background galaxies to define the transformation 
from each exposure into the reference frame. Although 
all stars can be fit with the same PSF, each galaxy has its 
own unique distribution of light and must be measured 
with its own template. An additional complication is 
that some images were taken at different orientations or 
with different instruments, such that the PSF for a given 
object might be different from epoch to epoch, affecting 
the observed distribution of light. For this reason, we 
construct templates for each galaxy with a deconvolved 
model, such that when we fit the template to the pixels in 
an image, we first convolve it with the PSF appropriate 
for that location in that detector. This approach ensures 
that the stars (which are simply delta functions when de¬ 
convolved) and galaxies are measured consistently within 
an exposure. 

We construct the star and galaxy templates from the 
exposures in the first-epoch data set. The transforma¬ 
tions (based on the stars) enable us to accurately map 
the pixels of all exposures into the reference frame. For 
each source, we collect all the pixels that map to within 
the 11x11 region about each source (galaxies and stars). 
We then use iterative forward modeling to construct a 
deconvolved template for each source. The template 
for the stars is simply a delta function. The template 
for the galaxies is a smooth empirical image that, when 
convolved with the PSF, best described the distribution 
of light in the contributing exposures. We allowed the 
galaxies to have an additional point source at their cen¬ 
ters, but most did not need one. 

We then use these templates, convolved with the ap¬ 
propriate PSF, to measure consistent positions for every 
source in every exposure in each epoch. We examine the 
fitting residuals for each source and reject the measure¬ 
ments that were clearly contaminated by cosmic rays or 
unidentified warm pixels. Measured positions are in the 
raw frame and must be corrected for geometric distor¬ 
tion, as mentioned above. Next, we use the positions for 
the galaxies to define the linear transformation from the 
distortion-corrected frame of each exposure into that of 
the first exposure of the first epoch. We examine the 
transformation residuals for each galaxy and reject those 
galaxies with unreliable and inconsistent positions. This 
leaves us with 10 to 15 high quality galaxies for each field. 
From the residuals, we infer that the transformations are 
typically good to ±0.93 mas, which corresponds to a 0.31 













Trajectories for 16 Runaways and HVSs 


5 


TABLE 1 

OBSERVED AND DERIVED HVS PROPERTIES 


ID 

RA, Dec 
(J2000) 

'^helio 

(km s“^) 

TeS 

(K) 

logs 

(cm s“^) 

(km s“^) 

Mg 

(mag) 

90 

(mag) 

^helio 

(kpc) 

Mra, / 
(mas y 

HVS 1 

9:07:44.99, +02:45:06.9 

831.1 ± 5.7 

11192 + 450 

3.91 + 0.20 

158 + 36 

-0.36 + 0.31 

19.69 + 0.023 

102 + 15 

+0.08 + 0.26, 

HVS 2 

9:33:20.87, +44:17:05.5 

708.0 + 15.0 

44561 + 675 

5.23 + 0.12 


+2.22 + 0.30 

18.56 + 0.013 

19 + 2.6 

-7.33 + 0.58, 

HVS 3 

4:38:12.77, -54:33:11.9 

723.0 + 3.0 

23000 + 1000 

3.95 + 0.10 

55 + 2 

-2.57 + 0.30 

16.36 + 0.20 

61 + 10 

+0.52 + 0.58, 

HVS 4 

9:13:01.00, +30:51:19.9 

600.9 + 6.2 

14547 + 598 

4.15 + 0.21 

77 + 40 

-0.71 + 0.33 

18.34 + 0.023 

64 + 9.8 

-0.23 + 0.36, 

HVS 5 

9:17:59.47, +67:22:38.3 

545.5 + 4.3 

12000 + 350 

3.89 + 0.13 

133 + 7 

-0.67 + 0.25 

17.58 + 0.032 

45 + 5.2 

+0.55 + 0.61, 

HVS 6 

11:05:57.45, +09:34:39.4 

609.4 + 6.8 

12190 + 537 

4.30 + 0.23 

170 + 55 

+0.25 + 0.27 

18.94 + 0.024 

55 + 6.9 

+0.05 + 0.57, 

HVS 7 

11:33:12.13, +01:08:24.8 

526.9 + 3.0 

12000 + 500 

3.80 + 0.10 

55 + 2 

-0.95 + 0.26 

17.63 + 0.015 

52 + 6.4 

+1.00 + 0.82, 

HVS 8 

9:42:14.03, +20:03:22.0 

499.3 + 2.9 

11000 + 1000 

3.75 + 0.25 

320 + 60 

-0.69 + 0.40 

17.93 + 0.016 

53 + 9.8 

-0.82 + 1.16, 

HVS 9 

10:21:37.08, -00:52:34.7 

616.8 + 5.1 

11637 + 520 

3.84 + 0.21 

306 + 72 

-0.71 + 0.34 

18.64 + 0.023 

74 + 12 

-1.26 + 0.74, 

HVS 10 

12:03:37.85, +18:02:50.3 

467.9 + 5.6 

11278 + 524 

4.38 + 0.23 

37 + 60 

+0.65 + 0.24 

19.24 + 0.024 

52 + 5.8 

-1.07 + 0.36, 

HVS 12 

10:50:09.60, +03:15:50.6 

552.2 + 6.6 

12098 + 622 

4.62 + 0.28 

78 + 88 

+0.55 + 0.28 

19.63 + 0.024 

66 + 8.5 

-0.40 + 0.36, 

HVS 13 

10:52:48.31, -00:01:33.9 

569.3 + 6.1 

11054 + 775 

4.00 + 0.35 

238 + 43 

-0.10 + 0.40 

20.01 + 0.021 

105 + 19 

-0.90 + 0.38, 

B 434 

11:02:24.37, +02:50:02.8 

443.9 + 2.9 

10232 + 382 

3.93 + 0.20 

117 + 42 

+0.06 + 0.27 

18.00 + 0.016 

39 + 4.8 

-1.61 + 0.28, 

B 485 

10:10:18.82, +30:20:28.1 

408.1 + 4.8 

16167 + 542 

4.02 + 0.12 

88 + 69 

-1.36 + 0.30 

16.06 + 0.016 

30 + 4.3 

-1.66 + 0.52, 

B 711 

14:20:01.94, +12:44:04.7 

273.7 + 5.4 

11004 + 298 

4.27 + 0.15 

60 + 76 

+0.72 + 0.26 

16.92 + 0.016 

17 + 2.0 

-0.96 + 0.80, 

B 733 

14:49:55.58, +31:03:51.3 

348.8 + 2.3 

10522 + 301 

4.25 + 0.14 

240 + 28 

+0.87 + 0.24 

15.67 + 0.020 

9 + 1.0 

-1.77 + 0.79, 


Referen ces. — (1) [Bro wri_ej^_alj ll2005[L_f^ljHIrsch et^aP (1200511 .^^CAjEdehnajmetaTjjl^OOSIL (4) lBrowri_e1^j^ |200^); (5) lBrown_et_j ^ |[20CK^; (6 ) ^mwn et 
<120091) : (8) IPrzvbilla et alTjl2008al) : fOl lBonanos et^^j|200^ : (10) IPrz vbi^^^^j|2008R) : (11) f200^j 7 fl2) lBro^^^ al.l J2012al) : {!{ 


TABLE 2 

LIST OF HST OBSERVATIONS 


ID 

UT Date 

Instr 

Filter 

P.A. 

(deg) 

Exptime 

(s) 

HVS 1 

2006-10-11 

ACS 

F814W 

-80.2 

522 X 5 


2009-10-05 

ACS 


-82.9 

540 X 5 


2013-02-21 

WFC3 


-80.2 

612 X 8 

HVS 2 

2006-10-04 

ACS 

F814W 

-56.2 

532 X 4 


2009-11-03 

ACS 


-56.2 

560 X 4 


2012-11-22 

WFC3 


142.8 

627 X 8 

HVS 3 

2006-07-08 

ACS 

F850LP 

-145.2 

257 X 6 


2009-12-23 

ACS 


35.0 

300 X 6 


2012-09-01 

WFC3 


124.8 

388 X 12 

HVS 4 

2006-11-06 

ACS 

F814W 

-82.0 

390 X 5 


2009-11-07 

ACS 


-86.0 

410 X 5 


2013-02-13 

WFC3 


24.8 

618 X 8 

HVS 5 

2006-09-01 

ACS 

F814W 

-26.1 

199 X 7 


2009-08-30 

ACS 


-25.9 

200 X 7 

HVS 6 

2009-11-07 

WFC3 

F606W 

159.3 

290 X 6 


2012-11-19 

WFC3 


159.6 

466 X 5 

HVS 7 

2009-12-04 

WFC3 

F606W 

157.5 

475 X 3, 90 X 3 


2012-12-05 

WFC3 


157.8 

531 X 3,160 X 3 

HVS 8 

2009-11-17 

WFC3 

F606W 

139.4 

460 X 3,120 X 3 


2012-11-26 

WFC3 


138.8 

494 X 3, 200 X 3 

HVS 9 

2009-11-07 

WFC3 

F606W 

151.8 

290 X 6 


2012-11-25 

WFC3 


152.1 

366 X 6 

HVS 10 

2010-02-23 

WFC3 

F606W 

103.4 

379 X 6 


2013-03-06 

WFC3 


103.3 

609 X 4 

HVS 12 

2009-11-08 

WFC3 

F606W 

155.2 

627 X 4 


2012-11-30 

WFC3 


155.2 

613 X 4 

HVS 13 

2009-11-13 

WFC3 

F606W 

154.5 

619 X 4 


2012-11-29 

WFC3 


154.5 

613 X 4 

B 434 

2009-11-06 

WFC3 

F606W 

155.3 

580 X 3,120 X 3 


2012-11-22 

WFC3 


155.6 

536 X 3,160 X 3 

B 485 

2009-10-30 

WFC3 

F606W 

161.4 

563 X 3, 25 X 3 


2012-11-20 

WFC3 


147.8 

675 X 3, 35 X 3 

B 711 

2009-12-18 

WFC3 

F606W 

170.5 

526 X 3, 50 X 3 


2012-12-16 

WFC3 


170.8 

608 X 3, 75 X 3 

B 733 

2009-12-12 

WFC3 

F606W 

-173.6 

563 X 3,16 X 3 


2012-12-09 

WFC3 


-173.3 

653 X 3, 30 X 3 


mas yr“^ systematic uncertainty in the reference frame 
for a 3 yr time baseline. The systematic uncertainty is a 
noise floor to all of our measurements. 

Finally, we use these galaxy-based transformations to 
map the position for each star in each exposure into 
the reference frame, and then solve for the proper mo¬ 


tions. Mathematically, the proper motion calculation 
boils down to a linear least squares fit to the x and y 
positions vs. time of observation, which we know exactly. 
We convert the x and y pixel motions into mas yr“^ us¬ 
ing the appropriate camera pixel scale, image rotation 
angle, and time baseline of observations. Proper motion 
uncertainties are determined from the scatter in x and y 
positions at each epoch, and added in quadrature to the 
uncertainty of the galaxy reference frame. Total mea¬ 
surement uncertainties range from ±0.35 mas yr“^ for 
the case of HVS 1, which has a 6.36 yr time baseline of 
observations, to ±1.33 mas yr“^ for the case of HVS 7, 
which has only 8 useable reference galaxies. Our final 
proper motion measurements are in Table [T] 


3.3. Comparison with Previous Measurements 

Our proper motion measurements are a 6-fold im¬ 
provement over previous measurements. One of the 
most accurate proper motion catalogs available today 
is the UCAC4, w hich covers the entire sky to a depth 
of i? = 16 mag (|Zacharias et al.l [201311 . One of our 
stars, B 733, is bright enough to have a UCAC4 mea¬ 
surement. The UCAC4 total proper motion for B 733, 
MUCAC4,B 733 = 13.9 ± 6.6 mas yr“^, agrees to within 
1 .5(7 of our measurement Mr 733 = 4.1 ± 1.2 m as yr“^. 
Another catalog is PPMXL (|Roeser et al.l [201^ . which 
combines 2MASS and USNO-B astrometry. Only B 733 
is bright enough to have a 2MASS measurement, and the 
PPMXL total proper motion ^ppmxl.b 733 = 6.3 ± 6.8 
mas yr“^ is in perfect agreement with our measurement. 
The SDSS survey probes much deeper and, wheii com¬ 
bined with USNO-B astrometry (iMonet et ai.l2003ll . pro¬ 
vides proper motion measurements for all of our sample 
except HVS 3. The average SDSS-USNO-B proper mo¬ 
tion uncertainty for our stars is ±5.3 mas yr“^, and only 
one star, HVS 5, has a SDSS proper motion that differs 
from zero at greater than 1.5 (t significance. In contrast, 
our HST measurements have an average proper motion 
uncertainty of ±0.8 mas yr“^, and thus are 6.6 times 
more accurate. 
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Fig. 3. — Proper motions in context: assuming a fixed distance, radial velocity, right ascension, and declination for each star (those 
listed in Table [T}, we compare the measured proper motions (points with errorbars) to the locus of trajectories that pass within 8 kpc 
(green ellipse) and 20 kpc (blue ellipse) of the Galactic center. Trajectories that pass through the Galactic center are marked by a +. 
All calculations use the lKenvon et al.l 1120141 1 potential model. To evaluate the likelihood of origin requires that we account for all of the 
measurement errors, as seen in Figure HI 


3.4. Proper Motions in Context 

At first glance, the stars have small proper motions and 
thus largely radial trajectories. Twelve of the stars have 
proper motions statistically consistent with zero. Only 
HVS 2, B 434, B 485, and B 733 have proper motions 
that differ from zero at better than 3 (t significance. 

For distant stars, however, the reflex motion caused by 
the Sun orbiting the Milky Way can dominate the appar¬ 
ent proper motion. The direction and amplitude of solar 
reflex motion depends on the location and distance of the 
star. Thus, a correct interpretation of the proper motion 
measurements requires that we calculate trajectories. 

As a first step, we adopt a fixed distance, radial veloc¬ 
ity, right ascension, and declination (those listed in Table 
[T|) for each star, and then calculate each trajectory back¬ 
ward in time for all possible p roper motions t hat cr oss the 
Galactic plane. We use the iKenvon et al.l (I2014D three 


component bulge-disk-halo potential model for the tra¬ 
jectory calculations. This potential model uniquely fits 
observed mass measurements from the Galactic center to 
the outer halo, but was originally constructe d for a 220 
km s“^ circular velocity (jKenvon et al.l[2f)?)^ . Updating 
the disk mass Md = 6 x 10^° Mq and disk radial scale 
length Od = 2.75 kpc yields a flat rotation curve of 235 
km s“^ consist ent with the mos t recent circular velocity 
measurements (iReid et al.||2014D . We refer to this poten¬ 
tial model as the IKenvon et all (|2014D model, and use it 
in all of our trajectory calculations. 

Assuming a fixed distance, radial velocity, right ascen¬ 
sion, and declination (those listed in Table [IJ, Figure [3] 
compares our measured proper motions to the results of 
our trajectory calculations. The green and blue ellipses 
show the locus of proper motions with trajectories that 
pass within 8 and 20 kpc, respectively, of the Galactic 
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center. Statistically, all of our stars have proper mo¬ 
tions consistent with a Milky Way origin. Our strongest 
constraints are for the three closest stars that have sig¬ 
nificantly non-zero proper motions: HVS 2, B 711, and 
B 733. To evaluate the likelihood of a more exact ori¬ 
gin requires that we account for all of the measurement 
errors, in proper motion, radial velocity, and distance. 

4. ORIGIN OF HYPERVELOCITY STARS 

With proper motions in hand, we can now address 
the question of origin. We expect that the fastest un¬ 
bound stars, as discussed above, are most likely HVSs 
ejected by the MBH in the Galactic center. Lower veloc¬ 
ity stars can be explained by alternative scenarios, such 
as runaway ejections from the disk. Importantly, run¬ 
aways with unbound speeds are most likely to be those 
stars ejected in the direction of Galactic rotation from 
the outer disk dHeber et al.ll2008l iBromlev et al.l [2?)0^ 
iKenvon et al.l 120141 . Thus we can distinguish between 
the runaway and MBH ejection scenarios if we know 
where the trajectories of our stars cross the Galactic 
plane. 

The stars in our sample were selected on the basis of 
their extreme radial velocities. Their trajectories in the 
Galactic frame will therefore point away from the Sun 
unless they have proper motions that significantly devi¬ 
ate from the Solar reflex motion. 

In reality, the errors in proper motion are too large to 
pinpoint the location of ejection with the desired accu¬ 
racy. To put the measurement uncertainties in context, 
consider a star’s tangential velocity Vtan = 4.74d/r, where 
d is the heliocentric distance in kpc and ^ is the proper 
motion in mas yr“^. The typical star in our sample has 
a median distance of 50 kpc and a proper motion error 
of ±0.8 mas yr“^, thus an error in tangential velocity 
of ±190 km s“^. For the median total space velocity of 
650 km s“^, the uncertainty in the angle of trajectory 
is ±16°. Over a distance of 50 kpc, the uncertainty in 
position is thus ±15 kpc. 

While we would like to test whether an observed trajec¬ 
tory goes exactly through the Galactic center, the errors 
preclude locating a trajectory to any one point. Instead, 
we use the measurements to try to disprove a Galactic 
center origin. Our approach is to calculate the statistical 
consistency of the measurements with a Galactic center 
trajectory. We obtain the same answer to this question 
whether we calculate trajectories from the star backward 
in time or from the Galactic center forward in time, be¬ 
cause the proper motion and distance errors are the same 
either way. Measurement errors make either the Galactic 
center or the star look very blurry. 

We use a Monte Carlo calculation to account for all of 
the measurement uncertainties, and visualize the results 
in a distribution of Galactic plane-crossing locations. For 
each star, we draw 1,000,000 current velocities and dis¬ 
tances assuming that the measured proper motion, radial 
velocity, and distance have Gaussian random uncertain¬ 
ties. We then calculate each trajectory backward in time 
and record where it crosses the Galactic plane. 

Figure S] presents the resulting distribution of Galac¬ 
tic plane-crossing locations for our sample. A black star 
shows the present location of each star, and a solid black 
line shows its trajectory. Magenta and cyan ellipses show 
the 1(7 and 2a likelihood plane-crossing regions, respec¬ 


tively. We determine these regions by calculating the 
density of plane-crossing locations in bins of X and Y, 
and solving for the density level thresholds that contain 
68.26% and 95.44% of all the crossings within them. The 
contours are centered on the mode of the plane-crossing 
distributions, which can differ from the mean trajectory 
(black line). Drawing both a large distance and a large 
proper motion can result in a very large tangential mo¬ 
tion and thus an extreme plane-crossing location. We 
discard those trajectories that fail to cross the Galactic 
plane within the main sequence lifetime of each star. For 
reference, a 3 Mq star has a 365 Myr main sequence life¬ 
time in the Padova tracks; for the special cases HVS 2 
and HVS 3, we allow flight times of I Gyr and 100 Myr, 
respectively. These flight time constraints exclude fewer 
than 1% of trajectories for most of the sample. 

Because distance and proper motion are the most un¬ 
certain parameters, the size and shape of the Galactic 
plane-crossing likelihood regions are driven by the dis¬ 
tance and proper motion uncertainties. There is no vis¬ 
ible change in the Galactic plane-crossing likelihood re¬ 
gions if we use a different pote n tial rn odel, specifically 
the one based on iGnedin et al.l (l2005h and updated in 
iGnedin et al.l (|2014h . Similarly, a triaxial halo potential 
model, like that expected from cosmological simulations, 
changes the trajectories very little. 

Half of our stars have trajectories that formally cross 
the Galactic plane at i? > 20 kpc, outside of the Milky 
Way disk, but all of our stars include a part of the Milky 
Way within their la plane-crossing likelihood regions. 
Only one star, HVS 7, has a trajectory moving opposite 
to the direction of Galactic rotation. HVS 7 is also our 
most uncertain proper motion measurement. Because 
a counter-rotation trajectory requires an unphysically 
large Galactic plane ejection velocity, we consider the 
astrometric link between the short- and long-exposures 
of HVS 7 suspect. Many more stars, including HVS 1, 
HVS 4, HVS 10, and B 485, have very radial trajectories. 

We reject the Galactic center origin hypothesis if the 
Galactic center trajectory falls outside the 3a {pcc < 
0.0026) threshold of the distribution (Figure |4]). The 
value of PGC is determined as the fraction of all possible 
orbits with the density of plane-crossings below the value 
at the GC bin (X=0,Y=0). This respresents the consis¬ 
tency of the Galactic center trajectory with the measure¬ 
ment uncertainties. By this estimate of probability, 13 
of our stars are consistent with a Galactic center origin 
within about the 2a {pcc ± 0.046) confidence level. As a 
sanity check, we obtain essentially the same result if we 
ignore the potential model and simply compare a per¬ 
fectly radial trajectory against the observed tangential 
velocities. Tangential velocity error drives the proba¬ 
bility. Only the stars HVS 2, B 711, and B 733 have 
measurements inconsistent with a Galactic center origin 
at >3(7 {pGC < 0.0026) confidence: they are runaways. 
We discuss these objects in more detail below. 

5. INDIVIDUAL OBJECTS 
5.1. HVS 2 

HVS 2 is a helium rich sd O star with a 708± 15 km s“^ 
heliocentric radial velocity (iHirsch et al.ll200^ . The two 
proposed explanations fo r its extreme motion are an ejec- 
tion by the central MBH (|Hirsch et al.ll200^lPeretsll2009ll 
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Fig. 4.— Distribution of Galactic plane-crossing locations for each star, calculated assuming Gaussian random uncertainties in the 
measured proper motions, radial velocities, and distances. For context, the large blue circles show the R = S solar circle and R = 20 kpc 
outer edge of the disk. The Sun is at X = —8 kpc (small blue circle), the Galactic center is marked with +, and the disk rotates clockwise 
in this plot. Black stars show the present location of each star, and the solid black lines with arrows show their trajectories. Magenta and 
cyan ellipses are the 1- and 2-cr likelihood regions where the stars cross the Galactic plane. 


or a Type la super nova explosion (| Just ham et al .112001 : 

I Wang fc Ha^[2009ll . 

HVS 2 has a significant proper motion of (/iRA, MDec) = 
(—7.33 ± 0.58,2.28 ± 0.55) mas yr“^. Our measurement 
is based on three epochs of imaging and therefore has a 
very well understood error distribut ion. Given HVS 2’ s 
18.5 ± 2.6 kpc heliocentric distance (|Hirsch et al.ll200'^ . 
the observed proper motion corresponds to a 673 ±118 
km tangential motion. HVS 2’s tangential motion is 
quite similar to its radial velocity. Thus its space motion 
in the rest frame of the Milky Way is almost exactly 1,000 
km s“^(!). 

The direction of HVS 2’s proper motion takes its tra¬ 
jectory across the Galactic plane at i? = 76.5 kpc. This 
result implies that HVS 2 originates from the stellar 
halo, possibly l aunched by a Type la s uper nova explosion 
as proposed by lJustham et al.l (I2009D and iWang fc HanI 


(|2009ll . If HVS 2’s distance and proper motion are Icr 
smaller, however, its tangential motion is less extreme 
and it crosses the Galactic plane around R = 20 kpc (see 
Figure n. Thus an ejection from the disk is a viable 
possibility. 

A Galactic center origin for HVS 2, on the other hand, 
is ruled out at greater than 3a confidence. HVS 2’s 
closest approach to the Galactic center, given the mea¬ 
surement errors, is i? = 4.5 kpc among our 1,000,000 
Monte Carlo trajectory calculations. In other words, 
the extreme velocity of HVS 2 cannot be explained by a 
dynamical interaction via gravitational interaction with 
the central MBH. Instead, it is an example of a hyper- 
r unaway star. 

iWang fc HanI (|2009ll perform ejection calculations for 
helium star companions to white dwarfs that explode in 
Type la supernovae. In this model, the helium stars re- 
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ceive a kick perpendicular to their orbital velocity for 
total ejection velo cities of 500 - 650 km s“^. Similarly, 
iGeier et all (|2013ll calculate a 600 km s“^ velocity from 
their Type la supernova ejection model. Given these 
ejection velocities, the Type la supernova model cannot 
explain HVS 2’s 1,000 km s“^ motion with a disk ejec¬ 
tion; HVS 2’s trajectory is not even in th e direction of 
Galactic rotation (Figure |4]). iTaurisI (|2015ll ’s aysmmetric 
core-collapse supernova model can produce 1,000 km s“^ 
velocities for low mass stars. However, there would not 
be enough time for a low mass star to evolve into an sdO 
star before its core-collapse companion launches it out 
of the Galaxy. The remaining possibility is that HVS 2 
was launched by a Type la supernova from a halo bi¬ 
nary traveling at ^400 km s“^ in its current direction of 
motion. 

Alternativ e orig ins seem less likely. For example, 
lAbadi et al.l (12009(1 propose that the tidal disruption of 
a dwarf galaxy can explain unbound stars. This mech¬ 
anism requires the close peri-center passage of a fairly 
massive > 10^° Mq dwarf to produce unbound stars 
(jPiffl et al.ll2011ll . The non-Galactic-center trajectory of 
HVS 2, plus the absence of other unbound stars around 
it, would appear to rule out the dwarf galaxy tidal de¬ 
bris origin. Additional evidence for a Type la supernova 
explosion, perhaps found in the abundance pattern of 
HVS 2’s stellar atmosphere, would better support the 
supernova origin picture. 

5.2. B 711 and B 733 

B 711 and B 733 are both bound B-type stars at mod¬ 
est 9-17 kpc distances, and thus plausible candidates 
for being runaway B stars ejected from the disk. Alter¬ 
natively, they could be failed HVSs on bound trajectories 
from the Galactic center. Their radial velocities in the 
Galactic frame are 290 and 440 km s“^, respectively. In 
the absence of a positive usini measurement for B 711, 
it is also possible that B 711 could be an evolved low 
mass star (i.e. a hot horizontal branch star) orbiting in 
the stellar halo. The fact that B 733 is a rapidly rotating 
2.5 Mq main-sequence star with a 349 km heliocen¬ 
tric radial velocity, however, requires that it was ejected 
from a location in the Galaxy with recent star formation. 

We measure proper motions that point to a disk run¬ 
away origin for both stars. The stars formally cross the 
Galactic plane at i? = 10 kpc (Figure H]). A Galactic 
center origin is ruled out at greater than Sa confidence 
for both stars. 

The trajectories of B 711 and B 733 cross the Galac¬ 
tic plane at angles nearly perpendicular to Galactic 
rotation, however, implying disk ejection velocities of 
533 km s“^ and 441 km s“^. The supernova ejec¬ 
tion mechanism is able to achieve 500 km s“^ veloc¬ 
ities foi_mam^quence stars only in extreme scenar¬ 
ios ((Taurisl [201511 . The dynamical ejection mechanism 
can achieve 500 km s“^ velocities for main-sequence 
stars but requires 3-body inter actions with contact bina¬ 
ries containing 100 Mg) stars (IGvaramadze et al.l 1^0091 : 
IGvaramadze fc GualandrTsI 1201 ih . Stars with 100 Mq 

are rare and short-lived. Simulations suggest that per¬ 
haps 0.1 % of dynamical eject ions reach 500 km s“^ ve¬ 
locities (|Perets fc Subd[20ll) . The upshot is that we 
expect to find more Galactic center ejections than disk 
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Fig. 5. — a) HVS 3 position measurements plotted relative to the 
mean; we label epoch 1 red, epoch 2 green, and epoch 3 blue, b) 
Dec positions versus date of observation; solid line is the linear least 
squares fit to the data, c) Same as b) but for RA. d) Proper mo¬ 
tion (star with errorbars) compared to the locus of HVS 3 proper 
motions with trajectories that pass with 3 kpc of the LMC (red 
ellipse) and within 8 and 20 kpc of the Milky Way (blue ellipses). 
The ellipses are calculated for a fixed HVS 3 distance; accounting 
for the distance error, HVS 3 is consistent with both LMC and 
Milky Way origins at the Icr level. 

runaway ejections a t these speeds (iBromlev et al.l (20091 : 
iPerets fc Subdl2012l : iKenvon et al.ll2014l) . The fact that 
B 711 and B 733 are ~500 km disk runaway ejections 
is thus quite intriguing. 

Because extreme runaway ejections require massive 
stars with relatively short lifetimes, we expect B 711 and 
B 733 to have flight times similar to their stellar ages. 
Our proper motions correspond to trajectories with 20 - 
40 Myr flight times from the disk. Our spectroscopic log g 
measurements favor young ages for both B 711 and B 733 
(see Figured]), however the uncertainties are large and no 
statistically meaningful constraint is currently possible. 
High resolution echelle spectroscopy of B 711 and B 733 
would thus be very interesting. B 733 is clearly a main- 
sequence star on the basis of its rapid rotation, and its 
trajectory is thus evidence for an extreme stellar dynam- 
ical ejection from the disk like that seen for HD 27179 1 
(|Heber et al.ll2008h and HIP 60350 (lirrgang et al.ll201^ . 

5.3. HVSS 

HVS 3 is the unbound 9 Mp^ main-sequence B star 
near the LMC on the sky (|Edelmann et al.l l2005fl . If 
HVS 3 comes from the LMC, then its speed is ev- 
idence for a MBH hidden somew here in the LMC 
(jGualandris fc Portegies Zwartll20n^ . If HVS 3 comes 
from the Milky Way, then its speed and stell ar nature 
are evidence of a former binary HVS ejection (iLu et al.l 
(20071 : |Peretsll2009ll . Both origins are unlikely in terms of 
ejection rates. 

In 2010, we published a prop er motion for HVS 3 that 
pointed to a Milky Way origin (|Brown et al.ll20lO) . This 
measurement was based on two epochs of ACS imaging 
with a 3.46 yr time baseline. The readout direction of 
the ACS CCD was unfortunately aligned with the Milky 
Way-LMC direction, however, and added an additional 
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systematic uncertainty. We now have a third epoch of 
WFC3 imaging, obtained at an orientation angle 90 deg 
from the previous data sets, that doubles our time base¬ 
line. Our expectation is that HVS 3’s intrinsic motion 
now dominates the errors. We re-process and re-analyze 
the epoch 1 and 2 data in the same way as the epoch 3 
data, so that everything is in a common reference frame. 

Figure [5] presents the results. We plot in panel a) the 
position of HVS 3, relative to its mean position, as mea¬ 
sured in each individual image. Different epochs are iden¬ 
tified by color, and the scatter in positions reflects the 
underlying precision of our measurements. We observe 
that HVS 3 moves 10.6 ± 5.0 mas, or 0.27 WFC3 pixels, 
in 6.15 years. This motion is quantified in Figure panels 
b) and c), which plot the RA and Dec positions of HVS 3 
versus time. The solid lines in each panel show the linear 
least squares fit to the measurements. The proper mo¬ 
tion of HVS 3 is (^iRA, MDec) = (0.52 ± 0.58,1.65 ± 0.57) 
mas yr“^. 

As fate would have it, our HVS 3 proper motion cor¬ 
responds to a physical trajectory that passes in-between 
the Milky Way and LMC. Figure [S] panel d) is identical 
to Figure [3] except that we now draw the locus of HVS 3 
proper motions with trajectories that pass within 3 kpc 
of the LMC (red ellipse). A 3 kpc radius encompasses 
the full extent of t he LMC bar and all of the young clus- 
ters proposed by IGualandris fc Portegies Zwarti (|2007l l 
for the origin of HVS 3. We determine HVS 3’s distance 
to the LMC in Figure [5] assuming that the LMC has a 
mass of 2 X 10^° Mq moving o n an orbit that reprod uces 
the line-of-sight velocity from Ivan der M arell (1200111 and 
the proper motion from iKallivavalil et al.l ( 20061) in our 
Milky Way potential model. If we account for the uncer¬ 
tainties in HVS 3 and LMC distances, then both LMC 
and Milky Way origins are consistent with the measure¬ 
ments at the 1(7 level. 

The ambiguity on origin is driven by our epoch 3 mea¬ 
surements, which exhibit a larger shift, and scatter, in 
position than seen in epochs 1 and 2. This result is 
visually apparent in Figure [S] Fitting only the epoch 
1 and 2 measurements yields {fiRA, ^J■Dec) = (0.10,0.89) 
mas yr ~^, in Icr agreemen t with our previously published 
value (|Brown et al.ll20iO) . Fitting only epoch 2 and 3 
measurements, on the other hand, yields ^J-Dec) = 

(1.01,2.52) mas yr“^. Clearly, the third epoch drives the 
proper motion to larger values, and thus causes HVS 3’s 
trajectory to move away from the Milky Way in physical 
space. 

We repeat our measurement and analysis for HVS 3 
using only the 5 best galaxies and obtain essentially the 
same result. Our reference frame thus appears robust. 
We speculate that the F850LP filter, chosen to maximize 
S/N on the background galaxies while not saturating the 
bright blue star, may be partly to blame. The F850LP 
filter has a less well-calibrated distortion solution than 
the F814W and F606W hlters, which were used for all 
of our other observations. The other HVSs with three 
epochs of data all have less scatter between epochs: the 
proper motions derived from epochs 1 and 2 and derived 
from epochs 2 and 3 agree at the la level, except for 
the Dec motion of HVS 2 that differs by 2cr. Residual 
CTE systematic error is also a possible problem for the 
HVS 3 observations. Answering the question of HVS 3’s 


origin will ultimately require a proper motion measure¬ 
ment with a longer time baseline of observations, or else 
a better instrument like Gaia. 

6. CONCLUSION 

We present HST proper motion measurements for 16 
stars with extreme radial velocities, 12 of which are un¬ 
bound to the Milky Way. On the basis of spectroscopic 
stellar atmosphere fits, our sample consists of 15 main- 
sequence B stars and 1 helium-rich sdO star located at 
10 - 100 kpc distances in the stellar halo. We expect that 
the fastest stars are likely HVSs ejected by the central 
MBH, and that the lower velocity stars may be runaways 
ejected from the Milky Way disk. 

We process our images using the best geometric dis¬ 
tortion solutions, CTE corrections, and empirical point 
spread function fits. Our final proper motions have 
an average uncertainty of ±0.8 mas yr“^, a 6-fold im¬ 
provement over previous measurements. Twelve of our 
stars have proper motions consistent with zero, and thus 
largely radial trajectories. 

Given the uncertainties in proper motion and distance, 
the data allow for a wide range of origin locations. We 
calculate the statistical consistency of the data with a 
Galactic center trajectory in an attempt to disprove the 
Galactic center origin hypothesis. We find that a Galac¬ 
tic center trajectory remains consistent with the mea¬ 
surements for 13 of our stars within the 2a confidence 
level. Only the stars HVS 2, B 711, and B 733 are incon¬ 
sistent with a Galactic center origin at >3 ct confidence, 
and thus runaways. 

HVS 2 is an unbound sdO star whose trajectory points 
from the stellar halo, possibly expl ained by a Type la su ¬ 
per nova explosion as pro posed by lJustham et al.l (1200911 
and IWang &: HanI (1200911 . Its 1,000 km s“^ motion is 
in some tension with the supernova ejection model, but 
can be explained if HVS 2 was ejected from a halo orbit. 
B 711 and B 733, on the other hand, are B stars with tra¬ 
jectories that clearly point from the stellar disk. These 
two stars are thus runaway stars, and their trajectories 
provide strong evidence for ^500 km ejections from 
the disk. 

Our third epoch of imaging for HVS 3, the unbound 9 
Mq B star near the LMC, yields a larger prop er motion 
than previously measured (iBrown et al.l 1201011 and thus 
a trajectory further from the Milky Way. Accounting for 
the uncertainty in HVS 3’s distance, its trajectory is now 
equally consistent with a Milky Way and LMC origin. 

In the future, better constraints on HVS trajectories 
and origin will come from improved proper motion mea¬ 
surements. Doubling the amount of spectra will not sig¬ 
nificantly improve distance estimates, for example, but 
doubling the time baseline of imaging will in principle 
double the precision of the proper motions. Even bet¬ 
ter, in 2017 Gaia will begin releasing proper motions for 
all of the HVSs. Gaia's predicted precision ranges from 
±0.03 mas yr“^ for 16th mag stars like HVS 3 and B 733, 
to ±0.3 mas yr“^ for a 20th mag star like HVS 1, and 
with improved accuracy compared to our measurements. 
Gaia will thus provide interesting constraints for all of 
the known HVSs. 

If the unbound HVSs are indeed ejected from the 
Galactic center, we can use their trajectories to probe 
the shape and orientation of the Milky Way’s dark matter 
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halo (|Gnedin et al.ll200^ . HVSs are effective test parti¬ 
cles that traverse the Galaxy to ~100 kpc distances. If 
the Galactic potential is triaxial, as predicted by cold 
dark matter simulations, the present motion of HVSs 
must deviate from being precisely radial. With a suf¬ 
ficient number of HVSs in different directions on the sky, 
and proper motions accurate to better than 0.1 mas yr“^, 
it may be possible to measure the two axis ratios and 
three direction angles of the triaxial halo. 

Facilities: HST (ACS, WFC3); MMT (Blue Channel 
Spectrograph) 
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